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Abstract 

A flight test of a business jet aircraft with various ice accretions was performed to obtain data to 
validate flight simulation models developed through wind tunnel tests. Three types of ice accretions were 
tested: pre-activation roughness, runback shapes that form downstream of the thermal wing ice protection 
system, and a wing ice protection system failure shape. The high fidelity flight simulation models of this 
business jet aircraft were validated using a software tool called “Overdrive.” Through comparisons of 
flight-extracted aerodynamic forces and moments to simulation-predicted forces and moments, the 
simulation models were successfully validated. Only minor adjustments in the simulation database were 
required to obtain adequate match, signifying the process used to develop the simulation models was 
successful. The simulation models were implemented in the NASA Ice Contamination Effects Flight 
Training Device (ICEFTD) to enable company pilots to evaluate flight characteristics of the simulation 
models. By and large, the pilots confirmed good similarities in the flight characteristics when compared to 
the real airplane. However, pilots noted pitch up tendencies at stall with the flaps extended that were not 
representative of the airplane and identified some differences in pilot forces. The elevator hinge moment 
model and implementation of the control forces on the ICEFTD were identified as a driver in the pitch 
ups and control force issues, and will be an area for future work. 


Nomenclature 


AOA 

angle of attack 

AvSSP 

Aviation Safety and Security Program 

C L 

aircraft lift coefficient 

Cl, max 

aircraft maximum lift coefficient 

dF 

flap deflection 

ICEFTD 

Ice Contamination Effects Flight Training Device 

IPS 

ice protection system 

LWC 

liquid water content 

MAC 

mean aerodynamic chord 

OAT 

outside air temperature 

TLF 

thrust for level flight 
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V 

airspeed 

Vfe 

flap extension speed 

VMO 

maximum operating speed 

Vs 

stall speed 

a 

aircraft angle of attack 


I. Introduction 

In an effort to improve flight safety, NASA’s Aviation Safety and Security Program (AvSSP) set out 
to develop and validate complete aircraft test methods and flight simulation models with emphasis on 
pilot training and aircraft icing certification testing. Specific objectives were to 1) develop piloted flight 
simulations to evaluate potentially risky flying qualities prior to development and certification flight tests 
with ice shapes and 2) develop piloted flight simulation models for operational training through higher- 
fidelity flight training devices. To fulfill these objectives, NASA teamed with Cessna Aircraft Company 
and Bihrle Applied Research to develop and validate flight simulation models of a typical business jet that 
incorporates the effects of various forms of ice. Prior to this effort, NASA, Bihrle Applied Research and 
the Wichita State University developed and validated a flight simulation demonstrator for a turbo-prop 
commuter aircraft for pilot training (refs. 1 to 3). 

The general methodology to develop the icing flight simulation models required static and dynamic 
wind tunnel testing of a subscale, complete aircraft model with and without ice shapes to develop a database 
of force and moment coefficients over a wide range of angles of attack, sideslip and control surface 
deflections. For the business jet simulation, a 1/12-scale complete aircraft model was used to develop the 
database for a no-ice baseline and three types of ice accretions. These were pre-activation roughness, 
runback shapes that form downstream of the thermal wing ice protection system, and a wing ice protection 
system failure shape. Because of the significant reduction in both the geometry and the Reynolds number, a 
series of wind tunnel tests was conducted on wing panel models prior to the 1/12-scale complete aircraft 
model testing. These wing panel model tests examined the scaling relationship between the full-scale 
aircraft and the subscale model with the three icing cases (ref. 4). From this scaling study, the subscale ice 
shape size, position, and roughness were determined for the complete aircraft testing. The 1/12-scale 
complete aircraft model was then tested on a rotary force-balance at the Bihrle Applied Research Large 
Amplitude Multi-Purpose Tunnel to develop the icing effects data base (ref. 5). Using this database, high- 
fidelity flight simulation models for the baseline (no-ice) and three icing cases were developed using the 
Bihrle Applied Research simulation environment software called D-Six. A limited validation exercise was 
conducted on the baseline data using flight data provided by Cessna. 

Additional flight test data with ice shapes were needed to validate the simulation models with the 
various forms of icing. Thus, a flight test program was conducted in Wichita, Kansas during September 
2005 using a Cessna business jet to acquire the full-scale flight dynamics with the no-ice baseline and three 
icing cases. The intent of this report is to describe the flight test performed to gather full-scale icing effects 
data for simulation model validation, discuss and share results from the simulation validation process, and 
describe simulation implementation on a flight training device to note company pilot comments comparing 
actual flight characteristics to the simulation models as implemented on the NASA ICEFTD. 

II. Experimental Methods 

A. Test Aircraft 

The test aircraft was a low-wing, T-tail twin pod-mounted turbo-fan light business jet which was 
modified for Cessna’s development testing (fig. 1). The test aircraft was instrumented to measure 
accelerations, rates and attitudes, air data variables, control surface deflections, engine fan speeds and 
throttle positions, fuel quantities, and pilot forces. Data was recorded onboard and transmitted through 
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Figure 1 . — Business jet test aircraft. 

telemetry to a ground station for engineering test support personnel to monitor progress and provide 
guidance on the test maneuvers. 


B. Simulated Ice Accretions 

Three types of ice accretions were tested on the full-scale aircraft. These were 1) pre-activation 
roughness, 2) runback ice shapes that form downstream of wing thermal anti-ice system, and 3) wing ice 
protection system failure case. These are typical ice accretions considered in certification programs 
(ref. 6). The geometries of the pre-activation roughness and failure shape were obtained from LEWICE 
1.6 calculations using the conditions in table 1. The runback shapes were developed from flight test 
experience in natural and tanker icing encounters. All three shapes will be discussed in more detail below. 


TABLE 1.— ICING CONDITIONS FOR PRE -ACTIVATION ROUGHNESS AND WING IPS FAILURE 


Ice shape 

Aircraft, 
a (deg) 

V, 

KCAS 

LWC, 

g/m 3 

MVD, 

jam 

OAT, 

F 

Time, 

min 

Pre-activation 

roughness 

3.08 

160 

0.3 

15 

-4 

2 

Wing IPS failure 

3.08 

160 

0.6 

15 

14 

22.5 


1. Pre-activation roughness configuration 

Pre-activation roughness represented a brief, initial icing encounter in which the ice protection 
systems were not active. Therefore, all flight surfaces (wing, vertical and horizontal tails) had a thin layer 
of ice roughness. The general position and locations of the pre- activation roughness was determined using 
LEWICE 1.6 and the conditions in table 1. LEWICE provided the geometric position of the icing extent 
for the wing and tail surfaces. Because the pre-activation shape did not have appreciable ice thickness, it 
was simulated using 40-grit sandpaper. Figure 2 shows the test aircraft with 40-grit sandpaper (blue areas) 
applied with double-sided tape on the leading edges of the wing and tail surfaces. 

2. Runback ice shape configuration 

The runback ice configuration represented an ice buildup after encountering an icing condition, such as 
very high liquid water content (LWC), which overwhelmed the wing thermal anti-icing system. In this 
scenario, some of the impacting liquid water would run aft of the ice protection system (IPS) and freeze 
on the upper and lower wing surfaces, while the leading edge remained free of ice. The other flight 
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Figure 2. — Pre-activation roughness (40-grit sandpaper-blue) on wing and tails. 

surfaces were assumed to have standard ice accretions, i.e., 45-min hold ice shapes on unprotected 
surfaces and inter-cycle ice on the horizontal tail. 

With that scenario in mind, the aircraft configuration for runback ice consisted of an upper wing 
runback shape, a lower wing runback shape, 45-min LEWICE shapes on the unprotected parts of the wing 
and horizontal tail, and a 45-min LEWICE shape on the vertical tail. The pneumatic horizontal tail ice 
protection system was assumed to be operational, so 40-grit sandpaper was used on protected area of the 
horizontal tail to represent inter-cycle ice. The upper wing runback shape was represented by a solid 
hump approximately 4 in. chordwise and a maximum thickness of about 1/2 in. It was positioned on the 
upper wing just aft of the thermal ice protection system. It was made from a silicone-based material and 
mounted using adhesive and aluminum speed tape. Forty -grit roughness was applied to the shape to 
provide surface roughness. The lower wing runback shape was a forward-facing wedge shape 
approximately 2 in. chordwise and 1 1/2 in. in height. It was positioned on the lower wing just aft of the 
thermal ice protection system. It was made from foam and mounted using double-sided tape. Forty-grit 
roughness was applied as well. The ice shapes on the unprotected portions of the wing and tail were 
derived using LEWICE with the same conditions as the Failed IPS shown in table 1, but for a 45-min 
exposure. On this aircraft, the unprotected portions were the wing tips, the entire vertical stabilizer, the 
roots and tips of the horizontal stabilizer. All of the LEWICE-predicted ice shapes were fabricated from 
foam, mounted with aluminum speed tape and covered with 40-grit roughness. 

Figure 3(a) and (b) show the right wing upper surface and lower surfaces of the runback configuration 
respectively. Note that the large shape on the wing tip is the 45-min LEWICE predicted shape on the 
unprotected wing tip. 

3. Wing IPS failure ice shape configuration 

The wing IPS failure configuration represented the ice buildup on the wing after an ice protection 
system failure [typically this is the amount of ice accumulated in 22 1/2 min (one-half of the 45-min hold 
condition) (ref. 6)]. The other flight surfaces were assumed to have standard ice accretions, i.e., 45-min 
hold ice shapes on the unprotected surfaces and inter-cycle ice on the horizontal tail. 

Figure 4(a) and (b) show the right wing and horizontal tail for this configuration. The 22 1/2 -min 
LEWICE shape was mounted to the wing leading edge nearly full span. This LEWICE failure shape 
varied along the span due to changes in the wing geometry (chord length and twist) and local flow 
conditions. The large shape on the wing tip is the 45-min LEWICE shape on the unprotected part of the 
wing. In the background, the silver area on the vertical tail shows where the 45-min LEWICE shape was 
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Figure 3(a). — Upper wing runback shape. Figure 3(b). — Lower wing runback shape. 



Figure 4(a). — Wing ice protection system failure case. Figure 4(b). — Horizontal and vertical tail ice for the wing 


ice protection system failure case. 

mounted. Figure 4(b) shows the right horizontal tail with the 45 -min LEWICE shapes on the tip and root 
and the 40-grit sandpaper (blue) over the leading edge to represent inter-cycle ice. All of the LEWICE 
shapes were fabricated from foam and mounted using aluminum tape. Forty-grit roughness was applied to 
these ice shapes for surface roughness effects. 

C. Flight Test Maneuvers 

The test maneuvers were selected to acquire airplane response characteristics for validating the flight 
simulation models for the baseline and the three ice configurations described above. Particular interest 
was placed on the higher angle of attack (AO A) where icing effects were more apparent and on control 
doublets to verify control effectiveness. All maneuvers were conducted with a center of gravity near 
25 percent mean aerodynamic chord (MAC) and initiated from trimmed steady-level conditions. All 
maneuvers were performed for 0° flap, while a subset of the maneuvers were done at 15° and 35° flap. 
Table 2 in the appendix lists all the maneuvers and configurations tested. Below is a brief description of 
the maneuvers performed in this flight test. 

1. Idle-power stalls 

The aircraft was trimmed at 1.3 stall speed (Vs) for the specific flap and gear position required at idle 
power. Pilot recorded the trim event, and then applied longitudinal control to decelerate the aircraft at 
1 knot per second. Pilot provided an event marker at the start of stick shaker and stall buffet and 
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continued decelerating until either full aft column position was reached for 2 sec or aerodynamic stall 
break occurred. Stall recovery consisted of lowering the nose to reduce AOA and adding power. The 
purpose of this maneuver was to determine the aircraft C L curve as a function of angle of attack. 

2. Elevator, aileron and rudder doublets 

The aircraft was trimmed at 1.3, 1.5, and 1.7 Vs or flap extension speed (V F E)for the specific flap and 
gear position required at thrust for level flight. Pilot provided an event marker for the trim event, and then 
proceeded with a controls fixed elevator doublet. After the motion damped out, the pilot would check the 
trim, provide an event marker and perform an aileron doublet. After the motion damped out, the pilot 
would check the trim, provide an event marker and then perform a rudder doublet. The purpose of these 
maneuvers was to acquire data on control effectiveness and damping characteristics in longitudinal, 
lateral, and directional axes. 

3. Steady-heading sideslips 

The aircraft was trimmed at 1.3, 1 .5, and 1 .7 Vs or V F e for the specific flap and gear position required 
at thrust for level flight. After noting the trim event, the pilot would perform steady a heading sideslip by 
slowly applying rudder while cross controlling with aileron input to maintain heading. The maneuver 
would end when full rudder deflection, 180 lb pedal force, 15° of sideslip, or buffet was attained. The 
pilot would slowly decrease rudder deflection and then repeat the process in the opposite direction. The 
purpose of this maneuver was to obtain data for lateral/directional forces, moments, and control 
effectiveness. 

4. Roll rate 

The aircraft was trimmed at 1.3, 1.5, and 1.7 Vs or V FE for the specific flap and gear position required 
at thrust for level flight. After noting the trim event, the pilot established a steady turn to the left with a 
30° ha nk angle. Once the left turn was established, the pilot used 1/3 aileron deflection to the right to roll 
the airplane to 30° right ha nk angle. The steady right turn was established, and then the pilot input a 
2/3 aileron deflection to the left to roll the airplane to 30° left ha nk angle. Again the steady left turn was 
established, and then the pilot input a full aileron deflection to the right to roll the airplane to 30° right 
ha nk angle. The purpose of this maneuver was to obtain data for lateral damping and control 
effectiveness. 

5. Maneuvering characteristics 

The aircraft was trimmed at 1.3, 1.5, and 1.7 Vs or V FF for the specific flap and gear position required 
at thrust for level flight. After noting the trim event, the pilot would perform a wind up turn by steadily 
rolling and increasing back pressure on the column until 2 G or buffet was attained. After the wind up 
turn, the aircraft was returned to straight and level flight, and a series of 2 G pull-ups and push-over 
maneuvers to 0.5 and 0 G were conducted. Pilots commented on the controllability of the aircraft. The 
purpose of these maneuvers was to obtain load factor, elevator, and column force data during each test 
condition. 

6. Thrust transients 

The aircraft was trimmed at 1.3 Vs for the specific flap and gear position required at thrust for level 
flight. After noting the trim event, the pilot rapidly advanced the throttles to maximum continuous thrust. 
The pilot maintained altitude while the aircraft accelerated by pushing the column forward. The maneuver 
was completed when maximum operating speed (VMO) or V FE was reached. The purpose of this 
maneuver was to identify thrust effects on trim control positions and pilot forces. 

7. Flap extension 

The aircraft was trimmed at 1.3 Vs or V FF for the specific flap and gear position required at thrust for 
level flight. Pilot provided an event marker for trim and then extended the flaps while maintaining initial 
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trimmed airspeed by adjusting aircraft attitude. Pilot commented on the control forces needed to maintain 
airspeed as the flaps extended. The purpose of this maneuver was to identify flap effects on changes in 
trim control positions and pilot forces. 


D. Flight Tests 

The flight tests were conducted in a build up of increased performance and control degradation. The 
initial flight was the no-ice baseline and the full test matrix in table 2 was accomplished in one flight. The 
next configuration was the Pre-activation Roughness. Since this configuration had not been flown in any 
previous tests, a series of buildup flights was required. In the first pre-activation flight, the 40-grit 
sandpaper was attached to the horizontal and vertical stabilizers and the inboard 1/3 span of the wing. 
After successfully qualifying this intermediate configuration, 40-grit sandpaper was attached to the next 
1/3 span of the wing and was flight qualified. The last 1/3 span was then covered and the full test matrix 
was accomplished for the pre-activation configuration in 2 flights. The next configuration flown was the 
runback ice, which had been flown in a previous program, so no buildup was required. All sandpaper was 
removed except for the horizontal tail ice protected area. The 45-min LEWICE shapes were attached to 
the unprotected leading edges, and the runback shapes were attached to the upper and lower surfaces of 
the wing. The full test matrix was accomplished in one flight. The last configuration was the Wing IPS 
Failure case. The Runback shapes were removed and a buildup of the Failure ice shape in 1/3 span 
increments (similar to the Roughness) took place to qualify this configuration. On the third flight, the full 
test matrix was conducted for the Wing IPS Failure configuration. 

The entire flight test was conducted in 1 1 flights (including qualification flights) in a total of 
18.3 flight hours. From those flights, a total of 200 test maneuvers were performed. 

E. Flight Simulator Sessions 

Prior to the flight test program, the preliminary flight simulation models of this business jet were 
implemented on NASA’s ICEFTD. The ICEFTD is a portable flight simulator that was developed as a 
concept demonstrator for building icing effects into flight simulation models for pilot training (refs. 2 and 
7). The business jet models had only limited validation for the no-ice case and the three iced cases were 
essentially unverified. The ICEFTD was shipped to Wichita, and during the course of the flight test 
program, company pilots flew various maneuvers in the ICEFTD with the various ice formations 
(fig. 5(a)). Pilots provided verbal comments comparing the simulator to the flight. 



Figure 5(a). — ICEFTD session during flight test. 


Figure 5(b). — ICEFTD session after model update. 
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After the flight test program, the flight simulation models were updated and re-implemented in the 
ICEFTD. Company pilots and engineers flew a series of maneuvers in the ICEFTD with the updated 
simulation models and again provided verbal comments comparing the simulator to flight (fig. 5(b)). 
These comments will be discussed in the next section. 


III. Results and Discussion 

A. Simulation Model Validation 

Flight test data were parsed into discrete data files for each maneuver and subsequently analyzed 
using the “Overdrive” tool in the Bihrle Applied Research’s D-Six simulation software (ref. 8) to validate 
the simulation models. “Overdrive” enabled the validation of the simulation aerodynamic database 
against flight-extracted data using the process illustrated in figure 6. At each time slice, extraction of 
aerodynamic moment coefficients from the flight-recorded time history occurred as shown on the right 
side of figure 6. Angular rates were numerically differentiated to obtain the angular acceleration of the 
vehicle. After the removal of the inertial effects, the remainder was nondimensionalized to calculate the 
aerodynamic force and moment coefficients experienced during flight. Also, at each time step, flight- 
recorded states, such as angle of attack, angle of sideslip, control surface positions, etc., were used to 
drive the aerodynamic model. Each aerodynamic model element (i.e., pitching moment due to elevator, 
etc.) was stored and summed as prescribed in the aerodynamic model. By over-plotting the model 
predicted coefficients with the flight-extracted total coefficients, differences were easily identified. When 
discrepancies were found, the simulation database was refined. The data would be rechecked with 
Overdrive to verify the refinement. 


1. General flight dynamics validation 

The validation effort showed that, in general, the original business jet simulation models that were 
developed from the wind tunnel data correlated well with the flight records. The validation effort resulted 
in only minor adjustments to the simulation models. The general flight characteristics for low to mid-level 

a, p, Speed, Altitude, Linear 



Figure 6. — Overdrive block diagram. 
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angles of attack were exceptionally well correlated. To verify general flight characteristics, maneuvers 
such as the steady heading sideslip and control doublets were examined. Examples of the “Overdrive” 
results for the steady heading sideslip and elevator, aileron and rudder control doublets are shown in 
figures 7 and 8 respectively. 

Figure 7 shows the Overdrive results from a steady heading sideslip maneuver. The time histories in 
the top 3 plots show the flight inputs while the 4th plot shows the simulation-predicted (blue) and flight 
derived (red) side force coefficient (CY) for this maneuver. The 5th plot displays the components of the 
simulation-prediction of CY. Note the close match between the simulation results and the flight derived 
terms as the sideslip (betad) varies over the range of-15° to 15°. This is an example of a good match for 
lateral directional characteristics. 



Figure 7. — Steady heading sideslip overdrive result. 
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Figure 8. — Control doublets overdrive result. 


Figure 8 shows the Overdrive results from a series of control doublets. The first input was an elevator 
doublet, the second was an aileron doublet, and the third was a rudder doublet. The first time histories 
show the angle of attack and sideslip from flight. The remaining three plots show the pitching moment, 
rolling moment and yawing moment time histories. The red traces are the flight derived terms while the 
blue traces are the simulation predicted moments. Clearly the simulation models accurately represent the 
flight records for these maneuvers. This suggests that the simulation models accurately represent control 
effectiveness and long and short period flight dynamics. 
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2. Icing Effects Model Validation 

The primary objective of the flight test program was to acquire data for validating the icing effects 
simulation models for the various ice accretions. Prior to this flight program, these simulation models 
were unverified. Critical points to accurately simulate were the speeds and angles of attack where stall 
occurred and the aircraft motion throughout the stall and recovery. The icing effects are clearly seen in the 
lift and pitching moment coefficients at higher angles of attack where flow separation and ultimately stall 
occurred. Figure 9 shows the lift curves for the no-ice and 3 iced cases derived from idle-power stall 
maneuvers to provide an overall sense of the effect of these ice accretions on the aircraft performance. 

The No-Ice lift curve provides a baseline reference. The lift plateaus at C imax with no sharp break at 
the maximum angle of attack. With the Pre-activation Roughness, the Ci max and stall angle of attack are 
nearly the same as the No-Ice baseline, but has a rather sharp stall break, and no plateau as was observed 
in the baseline. The larger ice accretions (Runback and Failure cases) had a more substantial degradation 
on lift performance. In both cases, the stall angle of attack was reduced by about 3° and the C imax was 
reduced by about 0.2. The character of the lift curves also resembled the No-Ice baseline in the flat 
plateau region prior to maximum angle of attack. Similar changes in stall characteristics due to the ice 
accretions were measured for the flaps extended configurations. 

To validate the icing effect in the simulation models, the Overdrive tool was used on the stall 
maneuvers in the same way as in the steady heading sideslip and the control doublets. Figure 10 shows 
the Overdrive result for the No-Ice, flaps up idle-power stall maneuver. The second plot in this figure 
displays the simulation-predicted lift coefficient as a function of time as well as the flight-extracted lift 
coefficient. A good match is observed throughout the maneuver even at high angles of attack and at stall 
recovery. Likewise, Overdrive was used on the stall maneuvers for each ice configuration. Figure 1 1 
shows the results for the Pre -Activation Roughness, flaps up configuration. The second plot shows a good 
match between the simulation-predicted and flight-extracted lift coefficients throughout the maneuver 
including stall recovery. Lastly, figure 12 shows the results for the Failed IPS Ice, flaps up configuration. 
Except for a small offset (0.02) in Cl at the beginning of the maneuver, the simulation-predicted lift 
coefficient matched the flight-extracted lift coefficient very well, particularly at stall and what appears to 
be a secondary stall that occurred during the recovery. 



No-Ice 

Roughness 

Failure 

Runback 


Figure 9. — Flight derived lift curves for various ice accretions (dF = 0). 
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B. Simulation Model Implementation and Evaluation on the ICEFTD 

Using the Overdrive tool, all the stall maneuvers for each flap setting and ice configuration were 
similarly validated. When discrepancies were found between the simulation-predicted and the flight- 
extracted force and moment coefficients, the deficiencies in the model were identified and corrected. 
Generally, the models required only minor adjustments. The model coefficient that required the most 
work was the pitching moment near stall for the flaps down cases. The original simulation model tended 
to have a pitch up tendency at stall. This was corrected by extending the nose-down pitching moment 
curve to a slightly higher angle of attack for these cases. 

The business jet simulation models were implemented on the ICEFTD so that company test pilots 
could experience the change in flight characteristics prior to the flight with each ice shape. Two company 
pilots flew the ICEFTD with each ice accretion configuration and provided comments based on their 
experience with the aircraft. The comments were favorable on the general flight characteristics. The 
longitudinal short period response was well damped, and the roll response was a good match. The primary 
discrepancy was with a pitch up that occurred during stall maneuvers with the flaps deflected to 15° and 
35°. Another discrepancy was with column forces. The ICEFTD provided control feedback through the 
column using a stick-loader system. The column forces were modeled using elevator hinge moment data 
in D-Six which provided commands to the stick-loader based on the flight conditions and pilot input. The 
stick forces were adjusted during this session, and pilots felt it resolved the issue in general. However, the 
forces were not representative at stall. The pilots expected a rather large pull force required during stall, 
but experienced push forces required to maintain pitch attitude. The pilot forces at the stall point appeared 
to drive the pitch up that was observed. 
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Figure 10. — Idle-power stall overdrive result (no-ice, dF = 0°). 
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Figure 11. — Idle-power stall overdrive (roughness, dF = 0°). 



Figure 12. — Idle-power stall overdrive (failed IPS ice, dF = 0°) 
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After the flight test, the simulation models were validated using the Overdrive tool as described 
above. The updated simulation models were re-implemented onto the ICEFTD, and the company pilots 
and engineers were invited to evaluate them. A test matrix similar to the flight test was developed, and 
each pilot flew each iced configuration at three flap settings to evaluate the general flight and stall 
characteristics. In general, both pilots considered the simulation to be representative for the flaps-up 
cases. Longitudinal and lateral/directional responses were similar to the aircraft, and stall characteristics 
were accurately represented for the No-Ice and various iced configurations. With flaps deflected, 
however, the pilots identified unrealistic pitch up tendencies at stall. 

Through this evaluation session, the team identified the elevator hinge moment model as a potential 
driver for these pitch up tendencies. At the angles of attack near stall, the ICEFTD control column forces 
tended to lighten and reverse. If unchecked by the pilot, the column movement would drive the elevator in 
a nose up manner and increase the angle of attack to very high (greater than 30°) values. In flight, the pull 
forces in the column increased as the elevator deflected trailing edge up, so that at stall, the pilot only 
needed to reduce the pull force to allow the column to move forward. In the ICEFTD, the pilots needed to 
push forward on the column to break the stall before very high angles of attack were attained. The 
elevator hinge moment model used in the ICEFTD was developed from low-speed wind tunnel test data 
provided by the company. This wind tunnel test data accounted for changes in elevator hinge moment due 
to flap deflection, but did not account for changes due to ice accretions. Additionally, the validation 
process using Overdrive focused on the verifying the aerodynamic forces and moments, not the control 
hinge moments. Therefore the hinge moment model was essentially unverified during these evaluation 
sessions. As a result of these evaluations, an additional validation effort focusing on the elevator hinge 
moment is warranted. 

IV. Conclusion 

A business jet aircraft was flight tested to measure the aerodynamic effects that three types of ice 
accretion had on the flight characteristics. The flight test data were used to validate flight simulation 
models that were developed through wind tunnel tests of a sub-scale complete aircraft model with test 
conditions representing the same three ice accretions. The validation process utilized a tool called 
Overdrive that compared flight-extracted aerodynamic forces and moments to simulation-predicted forces 
and moments. If discrepancies were discovered in these comparisons, the simulation engineer isolated the 
source of the discrepancy and adjusted the simulation model database to correct the model. The 
simulation models were implemented on a fixed based flight training device called the ICEFTD to enable 
pilot evaluations of the simulation models. 

The flight simulation models were successfully validated using the data collected through this flight 
test program. Validation was demonstrated through the Overdrive process which showed good match 
between flight-extracted forces and moments and simulation-derived forces and moments for the baseline 
and all three iced configurations. The adjustments to the simulation model database to achieve the good 
match were considered minor. Therefore the process to develop the icing effects simulation database was 
considered successful. 

Implementing these simulation models for pilot evaluation on the ICEFTD revealed an essential 
component in the simulation that had not been verified initially. The pilot-in-the-loop evaluations 
identified pitch up tendencies at stall with the flaps extended. The likely driver for these pitch ups was the 
elevator hinge moment model. Elevator hinge moments were modeled and implemented in the ICEFTD, 
but the model and the implementation of the model in the ICEFTD had not been verified with the same 
rigor as the aerodynamic forces and moments. The significance to this finding is that although the 
aerodynamic models may be accurate, a pilot may regard the simulation as inaccurate because of the 
effect that control forces have on the simulation. As a result of this finding, an additional effort will be 
made to verify the hinge-moment model and control force implementation on the ICEFTD. 
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Appendix 


TABLE 2.— FLIGHT TEST MANEUVERS AND CONFIGURATION MATRIX 


Cond 

Procedure 

WT 

CG 

FL/GR 

| Speed 

Alt 

Thrust 

1 

Flap extension 

Opt 

25% 

0/up 

on 

10 to 15 k 

TLF 

2 

Flap extension 

Opt 

25% 

15/dn 

warn 

10 to 15 k 

TLF 

3 

Flap extension 

Opt 

25% 

0/up 

200 

10 to 15 k 

TLF 

4 

Flap extension 

Opt 

25% 

15/dn 

161 

10 to 15 k 

TLF 









5 

Roll rate 

Opt 

25% 

0/up 


10 to 15 k 

TLF 

6 

Roll rate 

Opt 

25% 

0/up 

EEm 

10 to 15 k 

TLF 

7 

Roll rate 

Opt 

25% 

0/up 

iSSS 

10 to 15 k 

TLF 

8 

Roll rate 

Opt 

25% 

15/up 

Oil 

10 to 15 k 

TLF 

9 

Roll rate 

Opt 

25% 

15/up 

EEQM 

10 to 15 k 

TLF 

10 

Roll rate 

Opt 

25% 

15/up 

1 v FE 

10 to 15 k 

TLF 

11 

Roll rate 

Opt 

25% 

35/dn 

IEEH 

10 to 15 k 

TLF 

12 

Roll rate 

Opt 

25% 

35/dn 

on 

10 to 15 k 

TLF 

13 

Roll rate 

Opt 

25% 

35/dn 

1 v FE 

10 to 15 k 

TLF 









14 

Thrust transients 

Opt 

25% 

0/up 


10 to 15 k 

TLF 

15 

Thrust transients 

Opt 

25% 

15/up 

ims 

10 to 15 k 

TLF 

16 

Thrust transients 

Opt 

25% 

35/dn 


10 to 15 k 

TLF 









17 

Maneuvering char 

Opt 

25% 

0/up 


10 to 15 k 

TLF 

18 

Maneuvering char 

Opt 

25% 

0/up 

— W 

10 to 15 k 

TLF 

19 

Maneuvering char 

Opt 

25% 

0/up 

OB 

10 to 15 k 

TLF 

20 

Maneuvering char 

Opt 

25% 

0/up 

225 

10 to 15 k 

TLF 

21 

Maneuvering char 

Opt 

25% 

0/up 

250 

10 to 15 k 

TLF 









22 

Pitch, roll, yaw doublets 

Opt 

25% 

0/up 


10 to 15 k 

TLF 

23 

Pitch, roll, yaw doublets 

Opt 

25% 

0/up 


10 to 15 k 

TLF 

24 

Pitch, roll, yaw doublets 

Opt 

25% 

0/up 

■EEB 

10 to 15 k 

TLF 

25 

Pitch, roll, yaw doublets 

Opt 

25% 

15/up 

OT¥ 

10 to 15 k 

TLF 

26 

Pitch, roll, yaw doublets 

Opt 

25% 

15/up 

KOI 

10 to 15 k 

TLF 

27 

Pitch, roll, yaw doublets 

Opt 

25% 

15/up 

1 v EE 

10 to 15 k 

TLF 

28 

Pitch, roll, yaw doublets 

Opt 

25% 

35/dn 

BIM 

10 to 15 k 

TLF 

29 

Pitch, roll, yaw doublets 

Opt 

25% 

35/dn 

021 

10 to 15 k 

TLF 

30 

Pitch, roll, yaw doublets 

Opt 

25% 

35/dn 

1 v EE 

10 to 15 k 

TLF 









31 

SHSS 

Opt 

25% 

0/up 


10 to 15 k 

TLF 

32 

SHSS 

Opt 

25% 

0/up 


10 to 15 k 

TLF 

33 

SHSS 

Opt 

25% 

0/up 

O.^ 

10 to 15 k 

TLF 

34 

SHSS 

Opt 

25% 

15/up 

OH 

10 to 15 k 

TLF 

35 

SHSS 

Opt 

25% 

15/up 

HASH 

10 to 15 k 

TLF 

36 

SHSS 

Opt 

25% 

15/up 

1 v FE 

10 to 15 k 

TLF 

37 

SHSS 

Opt 

25% 

35/dn 

OSH 

10 to 15 k 

TLF 

38 

SHSS 

Opt 

25% 

35/dn 

on 

10 to 15 k 

TLF 

39 

SHSS 

Opt 

25% 

35/dn 

1 v FE 

10 to 15 k 

TLF 









40 

Stall characteristics 

Opt 

25% 

0/up 


10 to 15 k 

Idle 

41 

Stall characteristics 

Opt 

25% 

15/up 

EEB 

10 to 15 k 

Idle 

42 

Stall characteristics 

Opt 

25% 

35/dn 

OSH 

10 to 15 k 

Idle 
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